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& ■ Abstract 

Constraints on the EoS for symmetric matter (equal neutron and proton numbers) at supra- 
saturation densities have been extracted from energetic collisions of heavy ions. Collisions of 
neutron-deficient and neutron-rich heavy ions now provide initial constraints on the EoS of neutron- 
rich matter at sub-saturation densities. Comparisons are made to other available constraints. 



> 
o 



o 



1 Introduction 



The nuclear EoS is a fundamental property of nuclear matter that describes the relationships between 
^ ■ the energy, pressure, temperature, density and isospin asymmetry 5 = (p n — p p ) / (p n + p p ) for a nuclear 
0> system. (Here, p n and p p are the neutron and proton densities, respectively.) It can be divided into a 
symmetric matter contribution that is independent of the isospin asymmetry and a symmetry energy 
term, proportional to the square of the asymmetry, that describes the entire dependence of the EoS on 
^ ■ asymmetry. Investigations that provide an improved understanding of this term will also provide an 
^ ■ improved understanding of masses [1, 2, 3], fission barriers, energies of isovector collective vibrations 
[1, 4], and the thickness of the neutron skins of neutron-rich nuclei [5, 6]. 

Macroscopic quantities of asymmetric nuclear matter exist over a wide range of densities in neutron 
stars and in type II supernovae [7, 8]. Experimental information about the EoS can help to provide 
improved predictions for neutron star observables such as stellar radii and moments of inertia, crustal 
vibration frequencies [7, 9], and neutron star cooling rates [3, 7] that are currently being investigated 
with ground-based and satellite observatories. For many of these observables, the absence of strong 
constraints on the symmetry energy term of the EoS engenders major theoretical uncertainties. Conse- 
quently, the goal of determining the EoS has been a major motivation for many X-ray observations of 
neutron stars [10]. Experimental investigations are needed to check whether constraints derived from 
future neutron star observations can be supported by laboratory measurements. 

The total energy per nucleon (i.e. the Equation of State (EoS)) of cold nuclear matter can be written 
as the sum of a symmetry energy term and the energy per nucleon of symmetric matter, 

E(p, 5) = E (p) + E 5 ; E 5 = S(p) ■ 5 2 , (1) 
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Figure 1: Constraints on the symmetric matter EoS from heavy ion collisions. The exper- 
imental constraints (cross-hatched and shaded regions) and the theoretical calculations are 
discussed in the text. 



where S(p) describes the density dependence of the symmetry energy term. Measurements of isoscalar 
collective vibrations, collective flow and kaon production in energetic nucleus-nucleus collisions have 
constrained the equation of state for symmetric matter, E (p), and the pressure, P = p 2 ■ (dE (p)/dp\ s ), 
at T = and densities ranging from saturation density to five times saturation density [11, 12, 13]. 

Fig. 1 illustrates some of these constraints graphically. The shaded region extending from p/p = 2 
to p/po = 4.5 shows the range of pressures in cold symmetric nuclear matter for EoS's that describe 
the available transverse and elliptical flow data [11]. The center of the shaded region extending from 
p/p = 1.2 to p/po = 2.2 shows the pressure - density relationship that best describes the available 
kaon production data [12]. The range of densities for the kaon constraint represents an educated guess 
based on available calculations [12]; systematic studies to establish this range of pressures have not been 
performed. The dashed line extending approximately from p/po = 1-2 to p/po = 1-7 extrapolates the 
pressures consistent with the Giant Monopole Resonance data to higher densities. Various theoretical 
curves labeled RMF:NL3 [14], Akmal [15], and Fermi Gas illustrate the larger range of theoretical pre- 
dictions for the EoS at supra-saturation densities. Experimental constraints can discriminate between 
predicted behaviors of the EoS at supra-saturation densities and show which are more likely. 

The extrapolation of the EoS to neutron-rich matter depends on S(p), which has comparatively few 
experimental constraints [5]. Many recent efforts to constrain the density dependence of the symmetry 
energy have focused on its behavior near saturation density. There, one may expand the symmetry 
energy, S(p), about the saturation density as follows, 

S(P) = S + + %(^) 2 + ... , (2) 

3 po 18 po 

where L and K sym are slope and curvature parameters at p . The slope parameter, L, is related to 
Po, the pressure from the symmetry energy for pure neutron matter at saturation density via L = 
3p \dS p /dp\ po = [3/p ]/po- The symmetry pressure, p , provides the baryonic contribution to the 
pressure in neutron stars at saturation density [3], where the energy of symmetric matter, E Q (p), 
contributes no pressure, and it is also related to the neutron skin thickness, 5R np , of neutron rich heavy 
nuclei including 208 Pb [5, 6]. 
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2 Experimental Probes of S(p) 



Both symmetry and total binding energy terms in the nuclear semi-empirical binding energy formulae 
reflect averages of E(p,S) over the densities of nuclei [1]. The values of surface and volume symmetry 
energy terms obtained over fits of such formula to measured masses provide some sensitivity to the 
density dependence of S(p) near saturation density [1, 2]. Giant resonances, low-lying electric dipole 
excitations and the difference between neutron and proton matter radii may also provide sensitivity to 
the density dependence of S(p) near saturation density [4, 6, 16, 17]. 

Nuclear collisions can transiently produce both sub-saturation and supra-saturation density vari- 
ations. The symmetry energy has been recently probed at sub-saturation densities in collisions via 
measurements of isospin diffusion [18, 19], and of double ratios involving neutron and proton energy 
spectra [20]. Future comparisons of the spectra and flows of neutrons vs. protons and negative vs. pos- 
itive pions can provide constraints on the symmetry energy at supra-saturation densities [21, 22]. Such 
observables reflect the transport of nucleons under the influence of nuclear mean fields and of the colli- 
sions induced by the residual interactions, both of which can be described by transport theory. Focusing 
on sub-saturation densities, Chen et al., obtained a reasonable description of isospin diffusion data with 
a symmetry energies of approximate form with =0.69-1.1 [21, 26] using a Boltzmann-Uehling-Uhlenbeck 
(BUU) transport model IBUU04. While the double ratios of neutron and proton energy spectra [20] 
can also provide constraints at sub-saturation densities, these data have not been reproduced by the 
IBUU04 model [21]. 



3 Constraints on S(p) 

Recently, we have explored how isospin diffusion and n-p double ratios depend on both So and L 
parameters of S(p) using the Improved Quantum Molecular Dynamics (ImQMD) transport model[27]. 
Three observables, the neutron/proton double ratio [20], the isospin transport ratio at beam rapidity 
obtained from the isoscaling parameters [18], and the isospin transport ratios at various rapidities 
obtained from 7 Li/ 7 Be mirror nuclei ratios [19] were compared to ImQMD calculations for different 
assumptions about the density dependence of the symmetry energy [27]. 

Detailed description of the model and its application to the double ratio data can be found in ref. 
[28]. For brevity, we limit our discussion here to the parameterization of the symmetry energy used in 
our calculations, which is of the form 

S(p) = %f) 2/3 + %fF, (3) 
2 po 2 po 

Unless indicated otherwise, the kinetic and potential parameters are C s ^ = 25MeV, C s>p = 35.2MeV 
and the symmetry energy at saturation density, So = S(po) = 30.1MeV. 

We first turn our attention to the interpretation of neutron/proton double ratio data within the 
ImQMD model [28]. This observable derives its sensitivity to the symmetry energy from the opposite 
sign of the symmetry potential for neutrons as compared to protons [21]. First experimental comparisons 
of neutron to proton spectra in ref. [20] used a double ratio in order to reduce sensitivity to uncertainties 
in the neutron detection efficiencies and sensitivity to relative uncertainties in energy calibrations of 
neutrons and protons. This double ratio, 

np/vn/v/ u n m/p (0 , dM n (l)/dE c . m . dM p (2)/dE c . m . 
DR(Y(n)/Y( P )) = R n/p (l)/R n/p (2) = dMp{l)/dE _ " d M n (2)/dE c . m . ' (4) 

is constructed from the ratios of energy spectra, dM n /dE c . m . dM/dEC.M, of neutrons and protons for 
two systems 1 = 124 Sn + 124 Sn and 2 = U2 Sn + 112 Sn with different isospin asymmetries. The stars in 



3 




Figure 2: Left panel: Comparison of experimental double neutron-proton ratios (star sym- 
bols), as a function of nucleon center of mass energy, to ImQMD calculations (lines) with 
different density dependencies of the symmetry energy, parameterized by 7« in Eq. (1). Right 
panel: A plot of x 2 as a function of ji. 



the left panel of Fig. 2 show the resulting neutron-proton double ratios measured at 70° < 0cm < H0° 
as a function of the center-of-mass (CM.) energy of nucleons emitted from central collisions [20]. 

We have performed calculations for these two systems. Within statistical uncertainties, the double 
ratio observable, DR(Y(n)/Y (p)), is nearly independent of impact parameter for 1 < b < 5fm. The 
lines in the left panel of Fig. 2 show double ratios, averaged over 6=1,2, 3/m, vs. the CM. energy of 
nucleons for 7, = 0.35, 0.5, 0.75, 1 and 2. The theoretical uncertainties in Fig. 2 are statistical. Despite 
the large experimental uncertainties for higher energy data, these comparisons rule out both very soft 
(ji = 0.35, dotted line with closed diamond points) and very stiff (7, = 2, dotted line with open diamond 
symbols) density-dependent symmetry terms. The right panel shows the 7, dependence of the total x 2 
computed from the difference between predicted and measured double ratios. We determine, to within 
a 2cj uncertainty, parameter values of 0.4 < 7^ < 1.05 corresponding to an increase in x 2 by 4 above its 
minimum near 7, = 0.7. 

Decreasing ji corresponds to higher symmetry energy at sub-saturation densities, which enhances 
the emission of neutrons. Similarly, decreasing 7^ corresponds to increasing DR(Y(n)/Y(p)). In the 
limit of very small 7, << 0.35, however, the system completely disintegrates and DR(Y(n)/Y(p)) 
decreases towards the free nucleon value N tota i / Z tota i = 1.2. As a consequence of these two competing 
effects, the double ratio attains maximum values at about 7$ ~ 0.7. 

The density dependence of the symmetry energy has also been probed in peripheral collisions be- 
tween two nuclei with different isospin asymmetries by examining the diffusion of neutrons and protons 
between them. Unless the two nuclei separate, this "isospin diffusion" would continue until the chemical 
potentials for neutrons and protons in both nuclei become equal. To isolate diffusion effects from pre- 
equilibrium emission, Coulomb effects and secondary decays, measurements of isospin diffusion compare 
"mixed" collisions, involving a neutron-rich nucleus A and a neutron-deficient nucleus B, to the "sym- 
metric" collisions involving A + A and B + B. The degree of isospin equilibration in such collisions can 
be quantified by rescaling the isospin observable X according to the isospin transport ratio Ri(X) [18] 

Rt(X) = 2 ■ X-iX A+A + X B+B ),2 (5) 

~ -&B+B 
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Figure 3: Left panel: Comparison of experimental isospin transport ratios obtained from 
the yield ratios of A = 7 isotopes (star symbols), as a function of the rapidity, ImQMD 
calculations (lines) for b = 6fm. Right panel: \ 2 analysis for h = Qfm (solid curve) and 
b = 7 fm (dashed curve) as a function of 



where X is the isospin observable. Adopting the order of the notation, A + B to designate A and B as 
the projectile and target, respectively, one expects Ri(A + B) = R^A + A) and Ri(B + A) = Ri(B + B) 
for R values extracted from emission near projectile rapidities in the absence of isospin diffusion. In 
the opposite extreme, Ri(A + B) = Ri(B + A) m if isospin equilibrium is achieved. 

Eq. 5 dictates that different observables, X, provide the same results if they are linearly related 
[19]. Experimental isospin transport ratios obtained with isoscaling parameters, a, of charged particles 
of Z = 3 — 8, Ri(a) and yield ratios of A = 7 mirror nuclei, Rj = Ri(X 7 = ln(Y( 7 Li)/Y( 7 Be))) are 
consistent, ie. Ri(a) = Rj, reflecting linear relationships between a, Xj, and the asymmetry S of the 
emitting source [19]. For emission at a specific rapidity y, we assume R(a) = R(5) = Rj to be valid, as 
has been confirmed experimentally [19] and theoretically for all statistical and dynamical calculations 
[23, 24, 25]. In the following, we calculate 5 from the asymmetry of the fragments and free nucleons 
emitted at the relevant rapidity, but our conclusions do not significantly change if fragments alone are 
used to calculate 5. For reference, we note that BUU calculations use the asymmetry of the projectile 
residues [18] to construct Ri(5). 

Experimental isospin diffusion transport ratios, Rf, plotted as stars as a function of rapidity in the 
left panel of Fig. 3, have been obtained using Eq. 5 and the yield ratios of 7 Li and 7 Be, [19]. It is 
estimated that the measurements include impact parameters ranging from 5.8 to 7.5 fm [19]. We have 
performed ImQMD calculations at impact parameters of b = 5, 6, 7, and 8fm. The lines in the left panel 
of Fig. 3 show ImQMD calculations at b = 6fm . The corresponding \ 2 analysis at b = 6 (solid curve) 
and 7 (dashed curve) fm in the right panel displays sharp minima. Using the x 2 criterion adopted 
previously, the analysis favors the region 0.45 < 7« < 0.95. A similar analysis of the experimental 
and calculated isospin diffusion transport ratios, Ri(a), has also been performed and favors the region 
0.40 < 7i < 1.0 see ref. [27]. 

These constraints on the exponent 7i depend on the symmetry energy at saturation density, Sq = 
S(po), as well as L and K sym . Within families of parameterizations for S(p), K sym is largely determined 
by L. The present observables are more sensitive to L than to K sym ; L can therefore be determined 
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Figure 4: Representation of the constraints on parameters So and L. The various lines 
represent constraints discussed in the text. Corresponding values of po are given in the right 
axis. 



more reliably than K sym . We have preformed a series of ImQMD calculations at b = 6/m with different 
values of 7j and So to locate the approximate boundaries in the So and L plane that satisfy the 2a 
criterion in the x 2 analysis of the isospin diffusion data. The two diagonal lines in Fig. 4 represents 
estimates in such effort. The sensitivity of constraints on So and L to differences between m* and m*, 
and to the in-medium cross- sect ions have not been fully explored; we note that BUU calculations for 
isospin diffusion show much more sensitivity to Sq and L that to these other quantities [29]. The dashed, 
dot-dashed and solid lines centered around So = 30.1MeV in Fig. 4 represents L values consistent with 
the analysis of DR(Y{n) /Y{p)) (Fig. 2), Ri(a) and R 7 (Fig. 3), respectively. The vertical line at 
So = 31.6MeV depicts the range of L values obtained in ref. [21, 26] from comparisons of IBUU04 
calculations to the measured isospin diffusion data for Ri(a). Constraints from the isoscaling analyses 
in ref. [30] are not included as we have some concerns about their self-consistency, see ref. [31]. 

4 Context and Summary 

We have included, in Fig. 4, other recent constraints in the density dependence of the symmetry energy. 
The lower box centered at So = 32MeV depicts the range of po values from analyses of Pygmy Dipole 
Resonance (PDR) data [4]. The upper box centered at So = 32.5MeV depicts the constraints reported 
in Ref. [2] from the analyses of nuclear surface symmetry energies. The Giant Monopole Resonance 
(GMR) data of ref. [16] does not directly provide L; however, ref. [21, 16] reports the GMR data to be 
consistent with the isospin diffusion analysis denoted by the line at So = 31.QMeV . The inset of Fig. 
4 shows the density dependence of the symmetry energy that results from combining our shaded region 
with the limiting values of So = 30.2 and 33.8MeV provided by the PDR data [4]. 

In summary, both isospin diffusion and double ratio data involving neutron and proton spectra 
have been consistently described by a QMD model. The analyses of all three observables provide 
consistent constraints on the density dependence of the symmetry energy. The results overlap with 
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recent constraints obtained pygmy dipole and mass data. Some shifts in the boundaries of the constraints 
can be expected with improvements in the precision of the experimental data and in the understanding 
of theory. Nevertheless, the observed consistency between the different probes suggests that increasingly 
stringent constraints on the symmetry energy at sub-saturation density can be expected as the data 
improve and the model dependencies of their interpretations become better understood. 
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